Digital  imaging  of  clear-sky  polarization 
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If  digital  images  of  clear  daytime  or  twilight  skies  are  acquired  through  a  linear  polarizing  filter,  they  can 
be  combined  to  produce  high-resolution  maps  of  skylight  polarization.  Here  polarization  P  and  normal¬ 
ized  Stokes  parameter  Q  are  measured  near  sunset  at  one  inland  and  two  coastal  sites.  Maps  that 
include  the  principal  plane  consistently  show  that  the  familiar  Arago  and  Babinet  neutral  points  are  part 
of  broader  areas  in  which  skylight  polarization  is  often  indistinguishably  different  from  zero.  A  simple 
multiple-scattering  model  helps  explain  some  of  these  polarization  patterns.  ©  1998  Optical  Society  of 
America 
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1.  Introduction 

Observations  of  clear-sky  polarization  have  a  vener¬ 
able  history  that  includes  such  19th  century  icons  of 
optics  as  David  Brewster  (1781-1868),  Jacques  Babi¬ 
net  (1793-1872),  and  Francois  Arago  (1786-1853).1 
In  his  pivotal  1871  paper  on  sky  color  and  polariza¬ 
tion,  Lord  Rayleigh  (1842-1919)  provided  a  succinct 
theoretical  explanation  for  skylight’s  maximum  po¬ 
larization  90°  from  the  Sun.  Rayleigh  did  not  men¬ 
tion  the  neutral  points  (points  of  zero  polarization) 
observed  in  the  Sun-zenith  plane  by  Arago,  Babinet, 
and  Brewster.  Although  Rayleigh  was  confident 
that  he  had  “disposed  of  the  polarization”  problem,2 
others  did  not  view  his  silence  on  the  neutral  points 
as  a  proscription.3 

In  fact,  the  extensive  literature  on  clear-sky  polar¬ 
ization  left  by  Rayleigh’s  contemporaries  has  been 
greatly  expanded  in  20th  century  studies  by  (among 
others)  Dorno,  Jensen,4  Neuberger,5  Sekera,6  Coul- 
son,7  and  Volz.8  Many  of  these  studies’  measure¬ 
ments  have  been  made  with  naked-eye  polarimeters 
either  held  by  hand  or  mounted  on  a  clinometer.8’9 
Yet  even  electronic  polarimeters  usually  yield  only 
spot  observations  of  polarization  averaged  over  their 
narrow  fields  of  view  (FOV).  This  restriction  has 
tended  to  limit  the  angular  extent  of  polarization 
studies,  which  often  are  confined  to  the  Sun-zenith 
(or  principal)  plane  or  to  the  zenith  itself.7’10 
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Thus,  in  the  past,  skylight  polarization  was  ana¬ 
lyzed  outside  the  principal  plane  only  occasional¬ 
ly- 11,12  Yet  now  we  can  routinely  measure 
polarization  almost  anywhere  in  the  sky  hemisphere 
by  using  digital  imaging  techniques — in  my  case,  dig¬ 
itized  color  slides.  Because  skylight  is  mostly  lin¬ 
early  (rather  than  circularly)  polarized,13-15  I  can  use 
a  dichroic  linear  polarizing  filter  as  an  analyzer. 
The  resulting  maps  of  polarized  sky  radiances  can  be 
(1)  acquired  as  quickly  as  photographs  can  be  taken, 
and  (2)  readily  manipulated  to  produce  maps  of  sky¬ 
light  polarization.  In  addition  to  their  flexibility, 
digital-image  polarization  data  have  unprecedented 
angular  detail  and  extent.  The  high  spatial  and 
temporal  resolution  of  these  maps  can  improve  sky¬ 
light  polarization  models  as  well  as  the  physical  in¬ 
sights  that  underlie  them. 

2.  Measuring  Skylight’s  Degree  of  Linear  Polarization 

Compared  with  such  seemingly  straightforward 
quantities  as  skylight’s  radiance  and  chromaticity, 
measures  of  its  polarization  may  appear  obscure  to 
newcomers.  Graphs  of  polarization  ratio  versus 
scattering  angle16  and  hemispheric  maps  of  Stokes 
parameters  and  of  (partial)  polarization  planes17  can 
at  first  seem  far  removed  from  simple  observations  of 
skylight  through  a  polarizer.  Of  course,  historically, 
skylight  polarization  was  first  quantified  by  such 
simple  means.  To  identify  skylight’s  neutral  points, 
observers  such  as  Brewster  and  Arago  needed  only  to 
interpret  visually  the  images  that  they  saw  through  a 
polarimeter  (e.g.,  the  description  of  the  Savart  polari- 
scope’s  use  in  Ref.  5). 18 

To  quantify  polarization  further,  we  now  use 
Stokes  parameters  to  define  various  measurable  de¬ 
grees  of  polarization.  For  Stokes  parameters  I,  Q, 
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U,  and  V,  the  total  degree  of  polarization  PT  is  given 

by19 

(Q2  +  U2  +  V2)1'2 

PT—  -  ,  (1) 


where  I  is  the  scattered  irradiance  and  Q,  U,  and  V 
are  irradiances  related  to  the  scattered  light’s  ellip- 
sometric  parameters.20  By  choosing  a  linear  polar¬ 
izer  as  my  analyzer,  I  can  measure  only  skylight’s 
degree  of  linear  polarization  P,  where 


P  = 


(Q2  +  U2)1'2 

I 


(2) 


thus  ignoring  skylight’s  minimal  degree  of  circular 
polarization,  V/I.21 

Although  we  can  use  digital-image  radiances 
rather  than  irradiances  in  Eq.  (2),  this  still  leaves  the 
practical  problem  of  defining  the  Stokes  parameters 
in  terms  of  those  radiances.  For  measurements  of  P, 
Fitch  et  al.  define  convenient  relationships  between  I, 
Q,  and  U  and  digital-image  radiances.22  In  their 
work,  they  made  four  digital  images  of  a  scene  with  a 
linear  polarizer  rotated  to  four  different  relative  po¬ 
sitions:  0°,  45°,  90°,  and  135°.  The  0°  direction  can 
be  arbitrary  if  we  do  not  need  to  know  the  azimuth  x 
of  the  vibration  ellipse  for  skylight’s  polarized  com¬ 
ponent.23  For  linearly  polarized  light,  the  vibration 
ellipses  are  in  fact  lines.  Note  that  ellipsometric  az¬ 
imuths  x  are  different  from  4>rel,  azimuth  relative  to 
the  Sun  (see  below).24 

At  corresponding  pixels  in  a  scene’s  four  digital 
images,  the  relative  radiances  L(0°),  L(45°),  L(90°), 
and  L(135°)  are  related  to  the  Stokes  parameters  by 

I  =  0.5[L(0°)  +  L(45°)  +  L(90°)  +  L(135°)], 

Q  =  L(0°)  -  L( 90°), 

U  =  L(45°)  -  L(135°).  (3) 

Results  from  Eqs.  (3)  are  used  in  Eq.  (2)  to  calculate 
linear  polarization  P.  Equations  (3)’s  four-image 
technique  is  indispensable  if  a  lens’  FOV  is  quite 
large,  such  as  for  a  fish-eye  lens.  Yet  for  the  much 
smaller  FOV  of  a  normal  35-mm  camera  lens  (—37°  X 
25°),  even  two  images  can  provide  useful,  although 
incomplete,  information  about  clear-sky  polarization. 
Furthermore,  the  completeness  of  Eqs.  (3)  comes  at 
the  cost  of  extra  image  processing  and  added  uncer¬ 
tainties  in  geometrically  registering  four,  rather  than 
two,  images.  As  a  final  caution,  note  that  skylight 
polarization  has  been  measured  with  any  number  of 
polarizer  directions,25  and  each  of  these  techniques 
will  produce  similar  (but  not  likely  identical)  P  val¬ 
ues.  Every  additional  polarizer  direction  increases 
the  information  available  for  calculating  P  at  each 
pixel  and,  in  principle,  improves  the  accuracy  of  P. 

What  can  be  learned  from  the  simpler,  although 
more  limited,  two-image  technique?  First,  two- 
image  measurements  have  historical  significance;  in 
the  past,  zero  Q  was  often  equated  with  neutral  P.4-26 
Thus  seeing  the  relationship  between  Q  and  P  pro¬ 


vides  insights  into  earlier  research.  Second,  Stokes 
parameters  Q  and  U  are  useful  by  themselves  in 
quantifying  skylight  polarization,  and  either  can  be 
determined  from  pairs  of  polarized  digital  images. 

Suppose  that  we  (1)  include  the  principal  plane27 
(4>rei  =  0°  or  180°)  in  the  FOV  of  a  camera  equipped 
with  a  normal  lens  and  (2)  rotate  a  linear  polarizer 
mounted  on  that  lens  so  that  its  transmission  axis  is 
horizontal.  Above  the  Arago  neutral  point,  x  is  also 
horizontal,  and  so  the  angle  between  x  and  the  trans¬ 
mission  axis  is  0°  there  (x  is  vertical  below  the  Arago 
point).  Call  radiances  measured  at  this  filter  orien¬ 
tation  L( 0°).  Rotating  the  filter  90°  makes  its  trans¬ 
mission  axis  vertical  and  yields  L(90°).  Given  that 
U  =  Q  tan(2y)  and,  in  principle,  x  forms  angles  of  0° 
or  90°  with  the  principal  plane,  then  U  =  0  for  all 
skylight  observations  there.28  Now  define  the 
Q-component  contribution  to  polarization  as  pQ, 
where 


L(0°)  -  L(90°) 

Pq  L( 0°)  +  L(90°)  ’ 

with  90°  and  0°  measured  relative  to  x-  Pq  can  also 
be  regarded  as  a  normalized  Stokes  parameter  Q, 
although  without  Q’s  additive  properties  for  different 
scattering  mechanisms.  Naturally,  Eq.  (4)  offers  the 
advantage  of  requiring  only  two  polarized  digital  im¬ 
ages  rather  than  Eqs.  (3)’s  four. 

Assuming  U  =  0  in  the  principal  plane,  then  pQ  is 
the  signed  equivalent  of  Eq.  (2)’s  P  [except  for  small 
differences  in  the  denominators  of  Eqs.  (2)  and  (4)]. 
At  the  single-scattering  level,  Eq.  (4)’s  directions  are 
measured  with  respect  to  the  scattering  plane.  For 
single  scattering  by  spheres,  Eq.  (4)  in  fact  defines  the 
signed  polarization.29  P  +  \pq\  wherever  U  is  non¬ 
zero,  although  P  and  pq  may  actually  differ  very 
little.  In  the  principal  plane,  pQ  and  P  are  nearly 
identical. 

Outside  the  principal  plane,  linear  polarization  P 
still  can  be  determined  unequivocally  from  two  filter 
directions  if  the  detector  has  a  very  narrow  FOV  (say, 
<scl°).  Rotate  the  polarizing  filter  to  obtain  the 
FOV’s  radiance  extremes  Lmax  and  Lmin  and  then 
calculate  P  as30 


^max  -*-/min  /  r\ 

P  =  -  -  -  (5) 

Amax  Lmin 

Note  that  Eqs.  (2)  and  (5)  both  determine  the  same 
linear  P,  although  they  use  different  measurement 
techniques.  For  linearly  polarized  skylight,  Lmax 
and  Lmin  will  occur  at  filter  directions  90°  apart. 
Once  again,  extrapolate  this  to  a  camera  equipped 
with  a  normal  FOV  lens  and  polarizing  filter.  Al¬ 
though  the  camera  lens  has  a  FOV  ;§>  1°,  individual 
pixels  within  that  FOV  receive  radiances  from  very 
small  solid  angles  (<1  X  10  6  sr).  Aim  the  camera 
at  some  arbitrary  point  in  the  clear  sky  and  rotate  the 
polarizer  to  obtain  the  brightest  image  as  measured 
by  the  camera’s  light  meter.  Call  this  the  Lmax  im- 
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Fig.  1.  Signed  polarization  as  a  function  of  scattering  angle  'F  for 
scattering  by  a  small  sphere  (approximately  molecular  scattering) 
and  by  an  ensemble  of  spherical  haze  droplets  whose  size  spectrum 
follows  the  Deirmendjian  haze-M  distribution  at  X  =  475  nm.  The 
haze  droplets’  real  and  imaginary  refractive  indices  are  1.5  and 
0.01,  respectively. 


age;  rotating  the  polarizer  by  90°  yields  the  Lmin  im¬ 
age. 

Many  pixels  in  the  Lmax  and  the  Lmin  images  ap¬ 
proximately  meet  Eq.  (5)’s  conditions  for  calculating 
P.  Because  there  are  only  two  fixed  polarizer  direc¬ 
tions  and  x  varies  across  the  lens’  FOV,  many  pixels 
would  require  a  slightly  different  pair  of  polarizer 
directions  for  calculating  their  true  P.  In  fact,  be¬ 
cause  UP  0  outside  the  principal  plane,  we  are  ac¬ 
tually  calculating  pQ  for  most  pixels.  Yet  wherever 
(1)  polarization  is  strong  and  (2)  Eq.  (4)’s  0°  direction 
corresponds  to  x,  then  pQ  —  P  because  conditions  (1) 
and  (2)  require  that  \U\  <sc  |Q|.  Furthermore,  I  have 
independently  confirmed  that  pq  —  P  when  P  is  large. 
I  did  this  by  comparing  maps  of  two-image  pQ  and 
four-image  P  in  the  same  clear  sky  at  90°  from  the 
Sun.  In  this  strongly  polarized  sky,  patterns  of  pq 
and  P  are  virtually  the  same  and  x  changes  less  than 
4°  across  my  normal  lens’  FOV.  However,  when  x 
changes,  so  will  Eq.  (4)’s  pQ  values,  an  issue  ad¬ 
dressed  below.  Not  surprisingly,  the  combination  of 
variable  x  and  two  fixed  polarizer  directions  can  gen¬ 
erate  negative  pQ,  even  outside  the  principal  plane. 

Some  changes  in  x  (and  thus  pq  s  sign)  can  be 
explained  by  haze  scattering.  Figure  l’s  top  curve 
shows  polarization  as  a  function  of  scattering  angle 
for  an  isotropic  molecule.  Here  polarization  is  al¬ 
ways  nonnegative,  and  Rayleigh’s  theory  predicts 
complete  polarization  at  =  90°.  In  fact,  molecular 
anisotropy  reduces  the  maximum  single-scattering 
molecular  P  to  —0.92- 0.94. 31  Figure  l’s  bottom 
curve  shows  the  usually  negative  polarization  that 
results  from  scattering  by  the  Deirmendjian  haze-M 
droplet-size  spectrum  at  wavelength  X.  =  475  nm.32 
When  polarization  by  the  haze  drop-size  distribution 
is  equivalent  to  that  by  some  mean-sized  droplet, 
then  Fig.  1  also  gives  haze’s  single-scattering  contri¬ 


bution  to  pq.  In  a  multiple-scattering  atmosphere, 
the  relative  contributions  of  molecules  and  haze  to 
the  observed  polarization  will  be  complicated  func¬ 
tions  of  their  respective  slant  optical  thicknesses  and 
4'.  Yet  wherever  haze  scattering  predominates  that 
by  molecules,  pq  will  be  negative. 

So  in  general  |/5q|  P  P,  especially  if  Eq.  (4)’s  0° 
direction  is  chosen  arbitrarily.  This  inequality 
makes  sense  because  the  two-image  pQ  [Eq.  (4)]  is 
based  on  half  as  much  information  as  the  four-image 
P  [Eq.  (2)].  Yet,  as  shown  above,  P  can  be  measured 
unambiguously  from  only  two  radiances  [Eq.  (5)], 
provided  that  the  detector’s  FOV  is  quite  small.  Al¬ 
though  a  camera  with  a  normal  lens  and  polarizer 
takes  in  a  larger  FOV,  individual  pixels  within  that 
FOV  still  approximately  meet  Eq.  (5)’s  criteria  of 
maximum  and  minimum  radiance.  Furthermore,  if 
we  (1)  align  Eq.  (4)’s  0°  direction  with  x  and  (2)  mea¬ 
sure  pq  in  the  principal  plane  or  in  regions  of  strong 
polarization,  then  P  and  pq  will  be  nearly  the  same 
because  U  is  negligibly  small. 

3.  Polarizer  Performance  and  Possible  Error  Sources 

A  prerequisite  for  measuring  P  and  Pq  is  accurate 
clear-sky  radiances.  My  sources  for  these  are  color 
slides  of  clear  skies  photographed  at  several  sites  and 
times  of  day.  For  each  scene,  I  take  multiple  slides 
with  a  35-mm  camera  that  is  mounted  on  a  tripod  and 
equipped  with  a  normal  lens  and  linear  polarizing 
filter.  The  filter  is  rotated  to  yield  the  various  po¬ 
larized  radiances  on  the  right-hand  sides  of  Eqs.  (3) 
and  (4).  The  slides  are  digitized  at  a  color  resolution 
of  24  bits  per  pixel,  and  algorithms  developed  earli¬ 
er33  are  used  to  calibrate  the  digital  images  colori- 
metrically  and  radiometrically.  These  images  yield 
relative  radiances  whose  angular  resolution  is  lim¬ 
ited  only  by  the  film;  a  resolution  of  1/65°  is  possible. 

Although  the  camera’s  FOV  is  fixed  in  each  scene, 
the  polarized  images  still  must  be  geometrically  reg¬ 
istered  after  digitizing.  After  registration,  residual 
errors  vary  slowly  across  the  digital  images,  but  typ¬ 
ically  they  are  1:600  pixels  or  —0.07°  for  the  two- 
image  technique  [Eq.  (4)]  and  approximately  twice 
that  for  the  four-image  technique  [Eqs.  (3)].  Then 
Pq  or  P  are  calculated  pixel  by  pixel  from  the  entire 
set  of  digitized  color  slides.34  At  each  pixel,  values 
are  either  plotted  as  a  gray-scale  pixel  or  stored  for 
later  quantitative  analysis.  Note  that  (1)  all  my  ra¬ 
diance  and  polarization  measurements  are  inte¬ 
grated  across  the  visible  spectrum  (400-700  nm),  and 
(2)  the  spectral  transmittance  T  of  my  polarizing  fil¬ 
ter  is  nearly  constant  ( T  —  0.32)  at  visible  wave¬ 
lengths.  Solar  glare  complicates  making  accurate 
polarization  measurements,  so  my  polarization  maps 
exclude  the  Sun.  I  begin  with  two-image  measure¬ 
ments  of  skylight’s  pQ. 

In  taking  the  photographs  that  produce  Fig.  2(a),  I 
oriented  the  polarizing  filter’s  transmission  axis  par¬ 
allel  [L(90°)j  and  perpendicular  [L(0°)j  to  the  vertical 
principal  plane.  In  Fig.  2(a),  this  plane’s  projection 
is  a  vertical  passing  through  <firel  =  180°.  Below  the 
Arago  point  at  <l>ret  =  180°,  the  principal  plane  and  x 
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Fig.  2.  Maps  of  clear-sky  pQ  at  the  coastal  site  of  Chesapeake  Beach,  Md.,  on  6  February  1996.  In  map  (a),  the  Sun  elevation  =  1.8°  at 
2220  UTC  and  azimuth  relative  to  the  Sun  4>rel  =  180°  at  map  center.  In  map  (b),  the  Sun  elevation  =  -0.46°  at  2234  UTC  and  the  Pq 
maximum  occurs  at  4>rel  =  270°.  Each  map’s  angular  size  is  —23.5°  X  36°,  and  every  map  gray  level  spans  A pQ  =  0.05;  pixels  with  |pe| 
<  0.002  are  colored  red. 


coincide.  For  skylight  largely  due  to  molecular  scat¬ 
tering,  x  rotates  at  most  by  —5°  at  a  given  view- 
elevation  angle  0„  in  Fig.  2(a).35  Thus  the 
magnitude  of  pq  at  cf>rel  other  than  180°  will  be  re¬ 
duced  simply  because  the  polarizer  is  slightly  mis¬ 
aligned  with  x-  By  how  much  does  this 
misalignment  reduce  pq? 

The  short  answer  is  “some  in  principle,  but  almost 
none  in  practice.”  To  quantify  this  answer,  I  first 
determine  my  polarizer’s  optical  performance.  In 
polarizer  parlance,  the  extinction  coefficient  H90  is 
the  combined  transmittance  of  initially  unpolarized 
light  by  two  linear  polarizers  whose  transmission 
axes  are  crossed  at  right  angles.  The  photographic 
polarizer  used  here  has  an  average  visible  wave¬ 
length  Hgo  =  6.341  X  10  4,  and  its  major  and  minor 
principal  transmittances  k1  and  k2  are  0.6423  and 
9.873  X  10~4,  respectively.36  Given  these  data, 
Mueller  matrix  calculations37  show  that  my  polarizer 
can  transform  unpolarized  light  into  linearly  polar¬ 
ized  light  whose  P  =  0.9969.  Call  this  polarization 
that  is  due  to  the  filter  itself  its  intrinsic  polarization. 
Because  of  the  clear-sky  maximum  of  P  ~  0.85,  my 
filter’s  performance  is  entirely  adequate  for  measur¬ 
ing  skylight  polarization.  The  same  polarizing  filter 


is  used  for  all  measurements,  including  those  of  Figs. 
2-4. 

Now  we  can  answer  the  question  posed  above: 
how  does  polarizer  alignment  affect  skylight  pq ?38 
Figure  5  shows  Mueller  matrix  calculations  based  on 
my  filter’s  intrinsic  polarization  and  on  the  angle  ii 
between  its  transmission  axis  and  x-  These  calcula¬ 
tions  yield  Eq.  (5)’s  degree  of  polarization  for  com¬ 
pletely  linearly  polarized  light  ( P  =  1)  observed 
through  my  polarizer  when  it  is  misaligned  0  degrees 
with  respect  to  x-  This  underestimation  of  the  light 
source’s  true  polarization  is  a  depolarization  that  is 
caused  both  by  the  polarizer’s  intrinsic  polarization 
(. P  =  0.9969)  and  by  its  rotation.  Call  this  underes¬ 
timation  d(-&). 

As  a  multiplicative  factor,  d(d)  can  be  used  to  cor¬ 
rect  polarization  underestimations,  provided  that  we 
know  x’s  orientation.  In  principle,  that  orientation 
is  either  vertical  or  horizontal  at  4>rei  =  180°  and  0°. 
In  Figs.  2  and  3  my  threshold  for  neutral  pQ  is  0.002, 
and  I  have  colored  pixels  red  in  which  \pq\  <  0.002. 
Now  suppose  that  L(0°)  were  measured  with  the  po¬ 
larizer’s  transmission  axis  rotated  5°  from  x-  Then 
the  observed  threshold  pq  is  actually  0.002/d(5°),  or 
0.002037 — a  negligibly  small  error.  In  other  words, 
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Fig.  3.  Maps  of  clear-sky  polarization pQ  at  Marion  Center,  Pa.,  on  31  August  1996.  This  inland  site  is  —425  km  from  the  Atlantic  Ocean, 
and  its  elevation  is  450  m  above  mean  sea  level.  In  map  (a),  the  Sun  elevation  =  —0.7°  at  2352  UTC  and  the  center  <f>rel  =  0°;  in  map 
(b),  the  Sun  elevation  =  8.7°  at  2259  UTC  and  the  center  4>rel  =  104°.  Each  map’s  angular  size  is  —24.3°  X  35.7°. 


because  (iff!)  is  a  multiplier,  its  greatest  arithmetic 
effects  are  on  large  pq  values,  and  at  most  these  will 
be  only  —1.9%  if  ft  <  5°.  Thus  d(ft)  does  not  signif¬ 
icantly  increase  or  decrease  the  number  of  near¬ 
neutral  red  pixels  in  Fig.  2(a).  Counterbalancing 
(and  perhaps  outweighing)  this  small  depolarization 
is  the  camera  optics’  self-polarization  due  to  refrac¬ 
tion.39 

Other  possible  error  sources  are  likely  to  be  self- 
correcting.  Consider  (1)  geometric  exposure  falloff 
on  the  film  plane  (e.g.,  vignetting)  and  (2)  the  fact 
that  rays  with  different  incidence  angles  follow  dif¬ 
ferently  slanted  optical  paths  through  the  filter’s  di- 
chroic  layer.  In  principle,  both  factors  produce 
radial  radiance  patterns  in  the  digital  images.  How¬ 
ever,  because  the  exposure  reductions  are  the  same 
fraction  of  the  film-plane  radiances  for  all  filter  ori¬ 
entations,  their  effects  in  Eqs.  (2)  and  (4)  should  be 
self-canceling.  In  any  event,  Figs.  2-4  do  not  show 
any  consistent  radial  patterns  of  polarization  or  de¬ 
polarization. 

4.  Two-Image  Analyses  of  Clear-Sky  Polarization: 

A  Coastal  Site 

In  Fig.  2(a),  we  look  ENE  over  Maryland’s  Chesa¬ 
peake  Bay  just  before  sunset  on  6  February  1996. 


This  figure  is  a  gray-scale  map  of  pq  as  measured  by 
Eq.  (4)’s  two-image  technique  (map  time  and  place 
are  further  specified  below).  At  Fig.  2(a)’s  center, 
4>rel  =  180°,  so  there  we  are  looking  in  the  principal 
plane  where  pq  =  P  (assuming  U  =  0  there).  Figure 
2(b)  shows  pq  a  few  minutes  later  at  the  same  site. 
Now  4>rei  —  270°  at  image  center,  and  the  compass 
direction  there  is  SSE.  Each  map’s  angular  dimen¬ 
sions  in  Fig.  2  are  —23.5°  X  36°.  Throughout  Figs. 
2-4,  each  gray  level  spans  A pQ  =  0.05,  and  represen¬ 
tative  pq  values  are  labeled  on  the  maps. 

As  expected,  polarization  is  near  zero  some  dis¬ 
tance  above  the  antisolar  horizon.  Figure  2(a)’s  red 
pixels  (|pq|  <  0.002)  vividly  illustrate  that  the  Arago 
point  has  many  neutral  and  near-neutral  neighbors 
both  in  and  immediately  outside  the  principal  plane. 
(Recall  that  pQ  is  only  the  Q-component  contribution 
to  total  polarization.)  In  fact,  neutral  bands  such  as 
that  in  Fig.  2(a)  have  been  observed  before.  Dorno’s 
1917  observations  of  skylight  polarization  at  Davos, 
Switzerland,  included  what  he  termed  “neutral 
lines,”  which,  properly  speaking,  are  only  lines  of  Q 
(or  pq)  =  0.4’40  As  Coulson  notes,  “the  term  neutral 
line  is  somewhat  misleading,  as  P  A  0  at  neutral  lines 
except  as  they  intersect  the  Sun’s  vertical  at  the  neu¬ 
tral  points,  where  indeed  P  =  0.”28 
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Fig.  4.  Comparison  of  Pq  [map  (a)]  and  P  [map  (b)]  for  the  antisolar  clear  sky  at  Annapolis,  Md.,  on  12  March  1997  at  2255  UTC.  In 
both  maps,  <j>rei  =  180°  is  marked  and  the  Sun  elevation  =  2.4°.  In  map  (a),  pixels  are  colored  red  for  \pq\  <  0.002;  in  map  (b)  the  criterion 
is  P  <  0.005.  Each  map’s  angular  size  is  —36.2°  X  23.8°. 


Figure  6  shows  a  meridional  profile  of  pq  (~=P)  in 
Fig.  2(a)’s  principal  plane.  An  Arago  point  occurs  at 
0„  =  11.7°  in  Fig.  6,  some  13.5°  above  the  antisolar 
point.  Pixels  are  averaged  2°  on  either  side  of  the 
principal  plane  in  calculating  Fig.  6’s  observed  pQ, 
and  its  error  bars  are  based  on  standard  deviations  of 
Fig.  7’s  radiances  within  the  4°  azimuthal  swath. 
Figure  6’s  low  Arago  distances  have  been  observed 
before  over  water,  although  usually  for  higher  Sun 
elevations.41  The  boundary  between  negative  and 
positive  polarization  at  the  horizon  is  both  easy  to  see 
and  understand  here.  Skylight  reflected  from  the 
Chesapeake  Bay’s  planar  surface  toward  the  camera 
will  have  positive  polarization  for  many  different  in¬ 
cidence  angles  (here  a  typical  water  P  ~  0.15). 42 

Figure  2(b)’s  p(i  maximum  of  ~0.57  occurs  near 


image  center  at  c))rel  =  270°.  If  scattering  were  ex¬ 
actly  symmetric  about  the  principal  plane,  this  map 
would  be  identical  to  one  for  4>rel  =  90°.  However, 
the  irregular  pQ  patterns  around  4>rel  =  180°  and  0°  in 
Figs.  2(a)  and  3(a)  show  that  such  polarization  sym¬ 
metry  is  not  guaranteed.  In  Fig.  2(b),  the  0°  filter 
direction  was  the  one  that  gave  the  brightest  image  as 
measured  by  the  camera’s  light  meter.  In  this  L(0°) 
image,  the  filter’s  transmission  axis  was  rotated  —22° 
clockwise  from  the  vertical,  as  Fig.  2(b)’s  tilted  axis  of 
maximum  pq  suggests.  That  my  choice  of  L(0°)  and 
L(90°)  produced  large  pq  also  suggests  that  pQ  ~  P 
here.  Further  evidence  that  pq  ~  P  in  this  strongly 
polarized  region  comes  from  Fig.  2(b)’s  resemblance  to 
Bullrich’s  maps  of  polarization  maxima  observed  well 
outside  the  principal  plane:  Fig.  2(b)’s  maximum  is  at 
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polarizer  misalignment  ($) 


Fig.  5.  Depolarization  factor  d(fi)  as  a  function  of  the  misalign¬ 
ment  d  of  my  polarizing  filter’s  transmission  axis  relative  to  x-  If 
a  light  source  were  completely  linearly  polarized,  these  d(fi)  would 
equal  its  P  [Eq.  (5)]  as  measured  for  polarizer  misalignments  of  0° 
£  fi  <  5°. 


0„  ~  14.6°. 12  Pq  decreases  near  Fig.  2(b)’s  sides,  as 
expected  when  4'  approaches  180°  and  0°.  In  the 
haze  layer  just  above  the  horizon,  pQ  is  occasionally 
negative,  perhaps  because  negatively  polarized  scat¬ 
tering  by  haze  droplets  dominates  there  (see  Fig.  1  for 
70°  <  'VF  <  170°).  Also  evident  near  Fig.  2(b)’s  horizon 
is  a  steep  p()  gradient,  a  feature  caused  by  the  rapid 
increase  of  optical  thickness  with  zenith  angle.  The 
enhanced  multiple  scattering  that  results  will  rapidly 
depolarize  near-horizon  skylight  even  though  pci  is 


Fig.  6.  Clear-sky  polarization  component  Pq  («*P  here)  at  Ches¬ 
apeake  Beach,  Md.,  on  6  February  1996  at  2220  UTC  for  <f>rel  = 
180°  [see  Fig.  2(a)].  Each  error  bar  spans  two  standard  deviations 
of Pq  at  the  given  0„  and  is  calculated  from  the  standard  deviations 
of  Fig.  7’s  radiances.  The  above-horizon  rms  difference  between 
observed  polarizations  and  those  predicted  by  a  double-scattering 
model  =  0.0359.  The  model  wavelength  is  475  nm,  and  its  mo¬ 
lecular  and  aerosol  normal  optical  depths  are  0.15  and  0.05,  re¬ 
spectively. 


0.35  0.45  0.55  0.65 

scaled  radiance 


Fig.  7.  Polarized  clear-sky  radiances  looking  ENE  over  the  Ches¬ 
apeake  Bay  (Chesapeake  Beach,  Md.)  on  6  February  1996  at  2220 
UTC.  The  Sun’s  elevation  =  1.8°,  and  the  measurement  meridi¬ 
an’s  azimuth  relative  to  it  (4>rel)  is  180°.  Radiances  were  averaged 
2°  on  either  side  of  the  principal  plane,  and  the  L( 0°)  were  mea¬ 
sured  with  the  polarizing  filter’s  transmission  axis  perpendicular 
to  that  plane  (i.e.,  oriented  horizontal).  Multiplying  the  scaled 
radiances  by  LBcale  approximates  the  absolute  radiances  at  the  film 
plane,  including  the  effects  of  the  filter’s  average  spectral  trans¬ 
mittance  T. 


large  at  higher  0U.  Earlier  research  had  hinted  at 
such  polarization  gradients  near  the  horizon,43  but  dig¬ 
ital  imaging  reveals  them  in  unprecedented  detail. 

Figure  8  shows  meridional  radiance  profiles  of 
L(0°)  and  L(90°)  averaged  over  a  4°-wide  swath  that 
spans  Fig.  2(b)’s  polarization  maximum.  Except  at 
small  0,;,  Fig.  8’s  L(0°)  and  L(90°)  differ  much  more 


Fig.  8.  Polarized  clear-sky  radiances  looking  SSE  at  Chesapeake 
Beach,  Md.,  on  6  February  1996  at  2234  UTC.  The  Sun  eleva¬ 
tion  =  —0.46°,  and  <f>rei  =  270°.  Radiances  were  averaged  2°  on 
either  side  of  <|>rei  =  270°,  and  the  L(0°)  were  measured  at  the  filter 
orientation  that  produced  the  brightest  image  of  this  scene. 
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polarization  component  p^ 


Fig.  9.  Clear-sky  pQ  (~P)  at  Chesapeake  Beach,  Md.,  on  6  Feb¬ 
ruary  1996  at  2234  UTC  for  <j>rei  =  270°  [see  Fig.  2(b)].  The 
above-horizon  rms  difference  between  observed  polarizations  and 
those  predicted  by  a  double-scattering  model  =  0.0309.  The 
model  wavelength  is  475  nm,  and  its  molecular  and  aerosol  normal 
optical  depths  are  0.15  and  0.05,  respectively. 


Fig.  10.  Polarized  clear-sky  radiances  in  the  sunset  sky  near 
Marion  Center,  Pa.,  on  31  August  1996  at  2352  UTC.  The  Sun 
elevation  =  -0.7°,  and  4>rel  =  0°.  Radiances  were  averaged  2.5° 
on  either  side  of  the  principal  plane,  and  the  L(0°)  were  measured 
with  the  polarizing  filter’s  transmission  axis  perpendicular  to  that 
plane  (i.e.,  oriented  horizontally). 


than  they  do  in  Fig.  7.  Naturally,  that  large  differ¬ 
ence  appears  as  consistently  large  values  in  Fig.  9’s 
meridional  profile  of  pq  (~P).  Figure  9  also  shows  a 
rapid  decrease  in  pQ  above  ()„  =  32°,  and  Fig.  2(b) 
confirms  that  this  trend  occurs  at  all  4>rel  in  the  im¬ 
age.  Bullrich  also  observed  such  sharply  defined  po¬ 
larization  maxima.12 

5.  Two-Image  Analyses  of  Clear-Sky  Polarization: 

An  Inland  Site 

Figure  3(a)  maps  pq  in  a  mostly  clear  sunset  sky  near 
Marion  Center,  Pa.,  on  31  August  1996.  Similar  to 
those  of  Fig.  2,  the  angular  dimensions  of  Figs.  3(a) 
and  3(b)  are  each  ~24.3°  X  35.7°.  In  Fig.  3(a)  opti¬ 
cally  thin  clouds  cover  the  sky  within  9.5°  of  the 
topographic  horizon  (land  is  colored  green),  becoming 
optically  thick  only  at  0„  <  2.2°.  The  most  negative 
Pq  are  concentrated  not  in  the  optically  thick  clouds 
at  the  horizon  but  in  the  thinner  clouds  above  (note 
that  the  darkest  band  is  slightly  above  the  horizon). 
Scattering  geometry  and  optical  thicknesses  are 
much  different  than  in  Fig.  2(b),  yet  once  again 
greater  optical  thickness  is  associated  with  smaller 
(absolute)  pQ.  In  fact,  a  small  optically  thick  cloud 
atop  Fig.  3(a)’s  most  negatively  polarized  region  has 
some  positive  polarization  (see  points  just  below  the 
asterisks).  Recall  that  map  pixels  are  colored  red 
when  \pq\  <  0.002.  Isolated  red  pixels  within  the 
small  cloud  show  just  how  great  its  range  of  pq  is,  as 
it  is  embedded  in  a  region  where  \pQ\  >  0.05. 

A  neutral  band  is  evident  at  Fig.  3(a)’s  top,  as  is 
that  band’s  asymmetry  about  the  principal  plane.  A 
Babinet  neutral  point  exists  at  <t>rel  =  0°  within  the 
neutral  band  indicated  by  the  red  dots,  yet  dismissing 
all  of  its  red  neighbors  as  spurious  seems  unwar¬ 


ranted.  Although  polarization  models  may  define 
neutral  points  in  hypothetical  clear  skies,  in  real 
clear  skies,  pq  observations  at  best  show  only  neutral 
and  near-neutral  areas.  As  Fig.  4(b)  indicates,  even 
a  four-image  polarization  map  does  not  reduce  the 
Arago  area  to  a  point,  nor  will  allowing  red  pixels  only 
when  pq  =  0.0.  For  all  practical  purposes,  measure¬ 
ments  can  never  tell  us  when  pq  =  0.0, 44  and  thus  we 
can  never  insist  that  we  have  found  the  Babinet  (or 
Arago)  neutral  point  for  a  given  clear  sky. 

Figure  10  compares  the  polarized  radiance  profiles 
at  Marion  Center,  Pa.,  for  4>rei  =  0°.  As  in  Fig.  8,  Fig. 
10’s  L(0°)  were  measured  with  the  polarizer  trans¬ 
mission  axis  horizontal.  Here  L(90°)  exceeds  L(0°) 
at  most  0„,  making  pQ  negative.  The  L(0°)  and  the 
L(90°)  curves  cross  at  0„  =  28°,  the  elevation  angle  of 
a  Babinet  neutral  point.  Because  Fig.  3(a)’s  opti¬ 
cally  thick  small  cloud  is  silhouetted  against  the 
bright  twilight  sky,  it  appears  as  a  local  minimum  in 
the  radiance  profiles  at  0„  =  5.2°.  Yet  the  strong 
small-scale  depolarization  that  this  cloud  exhibits  in 
Fig.  3(a)  is  not  evident  in  Fig.  11  because  pq  s  under¬ 
lying  L( 0°)  and  L(90°)  have  been  averaged  over  5°  of 
azimuth  (see  Fig.  10).  However,  Fig.  ll’s  polariza¬ 
tion  profile  does  show  both  a  Babinet  point  and  the 
reduced  magnitude  of  pq  in  the  optically  thick  clouds 
near  the  horizon. 

In  Fig.  3(b)  we  see  pq  at  the  same  site  almost  1  h 
earlier,  only  now  the  map  is  nearly  perpendicular  to 
the  Sun’s  direction  (at  <brei  =  104°  here,  97°  <  T  < 
104°).  As  expected,  pQ  is  both  large  and  positive  in 
this  part  of  the  sky,  increasing  toward  the  left-hand 
side  of  the  map  (i.e.,  toward  c)>reI  =  90°).  As  in  Fig. 
2(b),  pq  ~  P  here  because  the  maximum-brightness 
L(0°)  image  produced  consistently  large  pQ.  Thin 
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Fig.  11.  Clear-sky  pQ  («*P  here)  near  Marion  Center,  Pa.,  on  31 
August  1996  at  2352  UTC  for  4>rei  =  0°  [see  Fig.  3(a)],  Each  error 
bar  spans  two  standard  deviations  of  Pq  at  the  given  0„  and  is 
calculated  from  the  standard  deviations  of  Fig.  10’s  radiances. 


clouds  near  the  horizon  eliminate  the  strong  pq  gra¬ 
dients  evident  above  Fig.  2(b)’s  horizon.  Several 
small  cumulus  appear  near  the  horizon  and  the  right- 
hand  sides  of  Fig.  3(b)’s  original  photographs.  These 
clouds  usually  are  local  minima  of  pq,  and  I  have 
marked  several  of  them  with  asterisks. 

However,  the  two  clouds  labeled  A  and  B  in  Fig. 
3(b)  are  both  local  minima  and  maxima  in  the  clear- 
sky  pq  field.  All  visible  parts  of  these  two  clouds 
were  sunlit,  and  the  optical  thickness  of  the  air 
column  between  the  camera  and  them  could  not 
have  varied  significantly  over  their  small  angular 
extents.  This  should  rule  out  any  differential  air- 
light  polarization  between  the  clouds’  tops  and  bot¬ 
toms.  The  optical  thicknesses  of  the  clouds 
themselves  do  not  appear  very  different  in  the  L(0°) 
and  L(90°)  images.  Given  that  T'  <  104°,  we  are 
not  near  the  cloudbow’s  polarization  maximum. 
Thus  we  have  no  ready  explanation  for  the  unusual 
polarization  patterns  in  clouds  A  and  B,  a  puzzle 
complicated  by  the  fact  that  other  small  cumulus 
nearby  are  consistently  depolarized  compared  with 
the  surrounding  clear  sky. 

Note  in  Fig.  12  that  0„  of  cloud  B’s  local  radiance 
maximum  shifts  ~1°  between  the  L(0°)  and  the 
L(90°)  images  (cloud  A’s  vertical  shift  is  much  small¬ 
er).  In  Fig.  13  we  see  the  corresponding  fluctuation 
in  cloud  B’s  pq  values.  Although  Fig.  12  helps  ex¬ 
plain  cloud  B’s  polarization  pattern  mathematically, 
it  provides  no  physical  insight  into  the  problem. 
And  even  if  four-image  polarization  P  eliminated 
cloud  B’s  local  polarization  maximum,  that  would  not 
explain  why  it  exists  for  pQ.  As  indicated  by  Fig. 
3(b),  all  other  clouds  in  Fig.  13  (including  A  for  c|)rel  = 
104°)  appear  only  as  depolarized  areas  in  the  clear- 
sky  pq  field.  None  of  this  complexity  was  evident 
during  photography,  when  I  saw  the  small  cumulus 


Fig.  12.  Polarized  clear-sky  radiances  near  Marion  Center,  Pa., 
on  31  August  1996  at  2259  UTC.  The  Sun  elevation  =  8.7°,  and 
c))rei  =  104°.  Radiances  were  averaged  2.5°  on  either  side  of  4>rei  = 
104°,  and  the  L(0°)  were  measured  at  the  filter  orientation  that 
produced  the  brightest  image  of  this  scene.  Here  small  cumulus 
usually  are  brighter  than  their  surroundings,  with  the  notable 
exception  of  cloud  D’s  L(0°). 


merely  as  annoying  contaminants  of  the  clear  sky. 
In  retrospect  these  clouds  pose  Fig.  3(b)’s  most  inter¬ 
esting  problem. 

6.  Simple  Double-Scattering  Model  of  Polarization 

Numerous  models  of  clear-sky  polarization  are  ex¬ 
tant,  10-11.4s  so  one  adds  to  the  list  with  caution. 
However,  as  a  first-approximation  comparison  of  ob¬ 
servation  and  theory,  my  earlier  model  of  clear-sky 
radiances  has  been  extended  to  include  pQ  (and  P  in 


Fig.  13.  Clear-sky pQ  (~P)  near  Marion  Center,  Pa.,  on  31  August 
1996  at  2259  UTC  for  <f>rei  =  104°  [see  Fig.  3(b)].  Although  most 
small  cumulus  here  are  local  Pq  minima,  note  that  cloud  B  is  both 
a  local  Pq  minimum  and  maximum. 
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the  principal  plane).46  The  extension  consists  of  us¬ 
ing  phase  functions  for  haze  and  molecular  scattering 
that  assume  polarized  rather  than  unpolarized  inci¬ 
dent  light.47  One  then  tracks  the  integrals  of  direct 
and  diffuse  skylight  contributions  along  a  specified 
slant  optical  path,  with  the  final  values  of 
perpendicular-  and  parallel-polarized  radiances  at 
the  observer  determining  pq.  The  model  yields 
quite  plausible  patterns  of  pQ,  including  a  band  of 
maximum  polarization  at  =  90°  and  neutral  bands 
such  as  those  seen  in  Figs.  2  and  3.  Whatever  its 
sins  of  incompleteness,  this  model  has  the  virtue  of 
being  a  straightforward  double-scattering  approxi¬ 
mation  to  the  radiative  transfer  equation. 

However,  my  goal  here  is  not  to  scrutinize  the  de¬ 
tails  of  an  admittedly  preliminary  model.  Rather,  I 
want  to  see  whether  it  can  begin  to  simulate  the 
detailed  pQ  observations  discussed  above.  Because 
this  simple  double-scattering  model  cannot  account 
for  polarization  asymmetries  about  the  principal 
plane  or  those  due  to  cloud  scattering,  I  have  not 
considered  the  Marion  Center  pQ  (see  Fig.  3).  How¬ 
ever,  by  using  the  Deirmendjian  haze-M  aerosol-size 
distribution  and  varying  the  aerosol  normal  optical 
depth  ja,  the  model  has  yielded  good  fits  to  the  Ches¬ 
apeake  Beach pQ  profiles  (see  Figs.  6  and  9).  At  A.  = 
475  nm,  the  best  fit  is  for  a  molecular  normal  optical 
depth  of  0.15  and  Ta  =  0.05.  In  the  antisolar  case 
(Fig.  6),  the  model’s  predicted  Arago  neutral  point  is 
1.7°  lower  than  that  observed,  and  the  two  polariza¬ 
tion  profiles  have  a  mis  difference  of  0.0359  above  the 
horizon.  Near  the  horizon,  however,  the  model’s  pq 
approaches  -0.2,  approximately  twice  the  observed 
polarization. 

The  much  better  near-horizon  agreement  seen  in 
Fig.  9,  coupled  with  the  fact  that  specularly  reflected 
sunlight  does  not  exist  there,  suggests  that  the  mod¬ 
el’s  Lambertian  surface  is  a  likely  source  of  Fig.  6’s 
errors.  Figure  9’s  rms  difference  of  0.0309  is  slightly 
smaller  than  Fig.  6’s,  with  the  largest  differences 
occurring  near  the  image’s  top,  where  the  observed 
polarization  decreases  rapidly.  A  major  failing  of 
the  model  is  that,  for  the  given  scattering  geometry 
(i.e.,  solar  elevation  =  -0.46°),  it  cannot  produce  p(j 
maxima  off  the  principal  plane,  contrary  to  what  we 
see  in  Fig.  2(b).  Possible  solutions  here  include  us¬ 
ing  a  wider  range  of  aerosol-scattering  phase  func¬ 
tions  and  examining  the  resulting  Mueller  matrices’ 
effects. 

Naturally,  the  comparisons  of  model  and  observa¬ 
tions  above  are  illustrative  rather  than  definitive. 
Duplicating  (or  even  approximating)  the  complex  po¬ 
larization  patterns  seen  in  Figs.  2  and  3  will  require 
a  much  more  sophisticated  model.  I  am  now  devel¬ 
oping  a  Monte  Carlo  model  that  will  explicitly  calcu¬ 
late  Mueller  matrix  transformations  of  incidence 
Stokes  parameters  at  each  scattering.  As  with  all 
Monte  Carlo  models,  one  can  ultimately  incorporate 
into  it  complex  spatial  distributions  of  scatterers  and 
surface  reflectance.  The  new  model  will  simulate 
patterns  of  P  that  could  be  compared  with  fish-eye 
maps  of  clear-sky  polarization.  Yet  even  the  present 


simple  form  of  the  double-scattering  model  offers  a 
level  of  realism  comparable  with  that  found  in  other 
models  (e.g.,  Ref.  1,  pp.  215-217).  As  an  instructive 
point  of  departure,  it  has  considerable  merit. 

7.  Two-  and  Four-Image  Maps  of  Antisolar  Clear-Sky 
Polarization 

By  how  much  do  the  two-image  pq  and  the  four- 
image  P  differ  in  a  given  scene?  Figure  4  shows  the 
complicated  relationship  between  the  two  polariza¬ 
tion  parameters  in  one  antisolar  sky  at  sunset.  The 
two  maps  are  of  the  same  clear  sky  photographed  at 
the  coastal  site  of  Annapolis,  Md.,  on  12  March  1997, 
when  the  Sun  was  2.4°  above  the  horizon.  In  Fig. 
4(a)  we  see  the  familiar  pQ  neutral  band  colored  in 
red  (|pq|  <  0.002),  which  once  again  is  asymmetric 
about  the  principal  plane.  As  expected,  negative  pq 
lie  below  this  band  and  positive  pQ  above  it.  As 
points  of  reference,  three  distant  radio  towers  near 
cj>rel  =  180°  form  dotted  red  lines  of  neutral  pq  that 
extend  from  the  neutral  band  to  the  horizon. 

In  contrast,  Fig.  4(b)’s  four-image  map  of  P  reduces 
the  near-neutral  band  to  a  roughly  circular  area  that 
straddles  the  principal  plane.  If  we  overlay  the  two 
maps,  the  near-neutral  P  region  falls  mostly  within 
its  pq  counterpart.  Although  that  overlap  is  ex¬ 
pected,  what  is  not  is  that  the  neutral  band’s  lowest 
0„  occurs  at  a  different  azimuth  than  does  the  P  neu¬ 
tral  area.  Furthermore,  P  maps  of  other  antisolar 
skies  at  sunset  show  circular  Arago  areas  several 
degrees  off  the  principal  plane.  Thus  we  have 
digital-image  confirmation  of  Brewster’s  (and  others’) 
naked-eye  observations  of  supposedly  anomalous 
neutral  points.48  Finally,  all  of  Fig.  4(b)’s  P  contours 
are  slightly  irregular,  suggesting  that  asymmetries 
in  pq  maps  are  not  artifacts  of  that  technique. 

Figure  14  shows  P  and  \pq\  averaged  over  a  5°-wide 
azimuthal  swath  that  is  centered  on  4>rei  =  180°  in 
Fig.  4.  As  above,  P  and  |/Jq|  are  essentially  the  same 
in  the  principal  plane,  even  though  their  averages 
include  polarization  data  from  just  outside  it.  In 
fact,  the  largest  difference  between  P  and  |/5q|  occurs 
near  the  pq  Arago  point  at  0„  =  9.3°.  Whatever  the 
higher  P  values’  physical  explanation  is  there,  their 
mathematical  explanation  is  clear  enough.  As  ex¬ 
pected  in  the  principal  plane,  the  average  U  is  much 
smaller  than  the  average  Q  (—4.6  times  smaller). 
Both  Q  and  U  change  sign  near  0,,  =  9.3°,  although  U 
is  negative  only  briefly.  However,  the  zeros  of  Q  and 
U  occur  at  different  0,,,  meaning  that  their  sum  in  Eq. 
(2)  is  never  zero.  Thus,  although  the  four-image  P  is 
very  small  near  pq  s  Arago  point,  that  average  P  is 
never  identically  zero,  even  if  we  smooth  over  a  nar¬ 
rower  azimuthal  swath.  However,  when  we  con¬ 
sider  individual  pixels  in  Fig.  4(b),  rather  than 
azimuthal  averages,  then  neutral  and  near-neutral  P 
values  abound. 

Note  that  the  four-image  technique  produces  an 
inherently  smoother  map  than  the  two-image  tech¬ 
nique.  Some  of  this  smoothing  stems  from  the  fact 
that  four  pixels,  rather  than  two,  are  added  to  form 
Eq.  (2)’s  normalizing  term  I.  Other  smoothing  is 
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Fig.  14.  Clear-sky  polarization  P  and  its  component  \pQ\  at  An¬ 
napolis,  Md.,  on  12  March  1997  at  2255  UTC  for  4>rel  =  180°  (see 
Fig.  4).  Both  P  and  pQ  were  averaged  2.5°  on  either  side  of  the 
principal  plane.  Each  error  bar  spans  two  standard  deviations  of 
P  at  the  given  0„.  Topography  that  extends  slightly  above  the 
astronomical  horizon  produces  the  large  variations  in  P  and  |p^| 
below  0„  =  2°. 


due  simply  to  the  four-image  technique’s  larger  pixel 
registration  errors.  This  four-image  smoothing  and 
its  resulting  polarization  reductions  are  evident  in 
Fig.  4(b),  in  which  red  pixels  have  all  but  disappeared 
below  the  horizon  and  in  the  radio  towers.  Less 
near-neutral  polarization  appears  in  Fig.  4(b),  de¬ 
spite  the  fact  that  its  threshold  for  red  pixels  has  been 
increased  to  P  <  0.005.  For  24-bit  images,  this 
larger  P  threshold  represents  a  difference  of  only  one 
or  two  gray  levels  in  each  color  channel  of  the  original 
digital  images.  Thus  we  still  are  close  to  the  thresh¬ 
old  of  quantization  errors.  In  any  event,  reducing 
the  neutral  P  threshold  to  0.002  (or  even  0.0)  simply 
yields  fewer  red  pixels  in  the  sky,  not  just  one. 

Thus  even  our  highest-resolution  analysis  of  P  will 
not  replace  Fig.  4(b)’s  near-neutral  area  with  a  single 
Arago  neutral  point.  This  should  not  alarm  us  any 
more  than  the  knowledge  that  a  pressure  minimum 
on  a  weather  map  corresponds  to  an  area,  rather  than 
a  point,  of  lowest  pressure.  Although  it  is  conve¬ 
nient  to  think  of  neutral  points  (or  even  neutral 
lines),  we  should  not  be  surprised  that  the  real  atmo¬ 
sphere  is  more  complicated  than  our  mental  (and 
computer)  models  of  it.  Similarly,  we  should  not  as¬ 
sume  that  neutral  P  outside  the  principal  plane  is 
somehow  an  anomaly.  The  atmosphere  is  a  tem¬ 
plate  for  our  models,  not  the  reverse. 

8.  Conclusions 

Digital-image  analysis  can  transform  clear-sky  pho¬ 
tographs  into  a  rich  and  reliable  source  of  polariza¬ 
tion  data.  This  new  data  source  not  only  has  very 


high  angular  and  temporal  resolutions,  but  it  also 
lets  us  see  polarization  patterns  over  large  areas  (e.g., 
Figs.  2-4).  These  maps  provide  several  new  in¬ 
sights.  First,  as  expected,  the  polarization  parame¬ 
ters  P  and Pq  do  have  different  spatial  patterns  in  the 
clear  sky,  but  clearly  these  are  related  to  each  other. 
As  Dorno’s  extensive  measurements  of  Q  suggest,  pq 
is  a  useful  polarization  measure.  Yet,  because  it  is 
based  on  less  information  than  P,  pQ  requires  espe¬ 
cially  careful  interpretation.  Second,  regardless  of 
whether  we  use  P  or  pQ,  the  Arago  neutral  point  (and 
presumably  others)  has  many  neutral  and  near¬ 
neutral  neighbors.  The  term  neutral  point  is  per¬ 
fectly  acceptable  so  long  as  we  remember  that 
skylight  polarization  actually  will  be  indistinguish- 
ably  different  from  zero  in  an  area  around  our  chosen 
neutral  point.  Third,  multiple  scattering  will  create 
strong  pq  and  P  gradients  in  the  near-horizon  sky  at 
right  angles  from  a  low  Sun.  These  gradients  are 
diminished  if  skylight  at  higher  0„  is  also  dominated 
by  multiple  scattering  (i.e.,  if  it  is  weakly  polarized). 
With  hindsight,  this  observation  may  seem  almost 
self-evident,  yet  Fig.  2(b)  provides  the  first  unambig¬ 
uous  confirmation  of  the  gradients’  existence. 
Fourth,  maps  of  both  P  and  Pq  indicate  that  skylight 
polarization  (and  radiance)  probably  should  be  as¬ 
sumed  to  be  asymmetric,  rather  than  symmetric, 
about  the  principal  plane.  As  convenient  as  the  lat¬ 
ter  assumption  is,  it  should  not  automatically  con¬ 
strain  skylight  models. 

Of  course,  we  need  not  fret  about  whether  a  model 
can  duplicate  every  irregularity  observed  in  a  given 
clear-sky  polarization  map.  More  important  is  that 
skylight  models  now  have  a  new,  and  much  more 
detailed,  standard  for  verisimilitude,  one  derived 
from  digital  imaging  of  clear-sky  polarization. 
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